Cosmological solutions with massive gravitons 
in the bigravity theory 

Mikhail S. Vo\ko^ 

Laboratoire de Mathematiques et Physique Theorique CNRS-UMR 6083, 
Universite de Tours, Pare de Grandmont, 37200 Tours, FRANCE 

We present solutions describing liomogeneous and isotropic cosmologies in the mas- 
sive gravity theory with two dynamical metrics recently proposed in larXiv: 1 109 .3515] 
and claimed to be ghost free. These solutions can be spatially open, closed, or flat, 
and at early times they are sourced by the perfect fluid, while the graviton mass 
typically manifests itself at late times by giving rise to a cosmological term. In ad- 
dition, there are also exotic solutions, for which already at early times, when the 
matter density is high, the contribution of the graviton mass to the energy density 
is negative and large enough to screen that of the matter contribution. The total 
energy can then be negative, which may result in removing the initial singularity. 
For special parameter values there are also solutions for which the two metrics ef- 
fectively decouple and evolve independently of each other. In the limit where one of 
the gravitational coupling constant vanishes, such special solutions reduce to those 
found in larXiv: 1107.5504] within the theory where one of the metrics is flat. 
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I. INTRODUCTION 

The currently observed acceleration of our universe l| is the main motivation of attempts 
to try to modify the theory of gravity, for example by giving a tiny mass to the graviton. This 
can effectively give rise to a small cosmological term leading to the late time acceleration 
The theory of massive gravity is not unique (see for a review) and there exist a number 
of its models, typically parameterized by two metrics, one of which is physical, (7^1/, while 
the other one is a fiat reference metric The coupling between the two is determined by 
a scalar function of g^'^fau- 



There 



Such models generically contain the negative norm ghost state in the spectrum 
is, however, a particular model, we shall call it RGT model, that could be special, since it 
is the o,Jy one that is ghost-ftee ■„ the deeoupUug hm.t g In add.t.o,, th,s .node, was 
recently claimed to be ghost-free in the full theory |6|. Even though the possibility to 
have a completely ghost-free massive gravity is sometimes disputed [7|, the RGT model is 
certainly interesting. Its solutions describing spatially closed, open and fiat homogeneous 
and isotropic self- accelerating cosmologies were obtained in js], jo]. 

Quite recently, the generalizations of the the GRT model were proposed, first by promo ;- 



ing the reference metric fu,u to be non-fiat [lOi], and next promoting it to be dynamical 



and it was argued that such generalized models remain ghost-free. In what follows we shall 



study the cosmological solutions within the bimetric theory of 

We find homogeneous and isotropic cosmologies which can be spatially open, closed, or 
fiat. For the generic parameter values they can be of two main types. First of all, there are 
solutions for which the universe expansion is driven at early times by the ordinary matter, 
while the graviton mass manifests itself only at late times by giving rise to a cosmological 
term. In addition, there are also exotic solutions, for which already at early times, when 
the matter density is high, the contribution of the graviton mass to the energy density is 
large and negative and screens that of the matter contribution. The total energy can then 
be negative, which may result in removing the initial singularity. In addition, for special 
parameter values, we find solutions for which the two metrics evolve independently of each 
other and the physical metric shows the late time acceleration. In the limit where one of 
the gravitational coupling constant vanishes, we shall call it the RGT limit, they reduce to 
those of RGT theory found in jsl. It turns out that the generic solutions do not reduce in 
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the GRT limit to any of the GRT theory, because both metrics remain then curved. 

In what follows we present a rather detailed analysis of these solutions as well as their 
relation to the GRT limit. We address, in particular, the question (for some reason not 
discussed in the literature) of how to derive the Lagrangian field equations in the theory. 
This issue is in fact not as trivial as it may seem, since the interaction between the two 
metrics is parameterized by j'^j^ = y/g'^'^fau, whose direct variation is problematic, since the 
variation (57'^^, does not commute in the matrix sense with 'y'^^. We show how this problem 
can be handled withing the tetrad approach. 

The rest of the paper is organized as follows. In SecJTl] we describe the bimetric general- 
ization of the RGT theory, the tetrad approach, the variation procedure and the reduction 
to the spherically symmetric sector. Solution with the decoupled metrics arising for special 
parameter values are described in SecJIIIt while SecJIVI presents a detailed analysis of the 
generic case. Yet one more class of solutions, arising due to a different way to fulfill the con- 
servation condition, is described in SecJVl while the relation to the GRT limit is discussed in 
Sec JVIl We summarize our results in Sec JVIII and list in the Appendix the energy- momentum 
tensor components. 

II. THE BIMETRIC GENERALIZATION OF THE RGT MODEL 

The theory is defined on a four dimensional spacetime manifold Ai parameterized by 
coordinates and equipped with two metrics g^u{x) and f^u{x). Their kinetic terms are 
chosen to be of the standard Einstein-Hilbert form, with the corresponding couplings G and 
TjG. The dynamics is determined by the action 

' = I (i« + ".^v^<''-^/^V^rf'x + S.-„ (1) 

where R and TZ are the Ricci scalars for gf^j, and /^y, respectively, while describes ordi- 
nary matter (for example perfect fluid) which is supposed to interact only with g^y. The 
interaction between the two metrics is defined by 

with 

K = K - , (3) 



4 



where 7^^^ is defined by the relation 

1^1% = g'^Uu , (4) 

and g^^ is the inverse of g^^ij. Apart from the gravitational coupling G, the theory contains 
three parameters 77, C3, C4. In the limit where — )■ and f^y is flat it reduces to the RGT 
theory [sl. 

The field equations are obtained by varying the action with respect to the metrics. A 
difficulty arises at this point, since varying the constraint (jlj) gives 

+ 1^5^. = Sg^^'U + g'^^Uy (5) 

and it is not obvious how to resolve this relation with respect to 5'^^^. One could of course 
consider 7^^^ as independent variables and impose the constraint (j4]) within the Lagrange 
multiplier method. However, the Lagrange multiplier enters then the equations as an aux- 
iliary field, and it is unclear how to get rid of it. Therefore, we adopt a different strategy, 
motivated by the approach of [sl. Let us introduce two tetrads e^j^ and oj^ defined by 
the conditions 

g^'^ = r^^^e^^e^ , f,. = vabco^co^ , (6) 

where tjab = diag(l, —1, —1, —1) is the Minkowski metric. We shall also need their inverse 
and such that e^e^ = 6^ and e'^e^ = 6^, similarly for w^. Each of the tetrads and 
is defined up to 6 local Lorentz rotation, so that equations contain a 12-parameter 
gauge freedom. Let us now require that 

eA^B^, = e'^ujAtM , (7) 

where ua^i = t^ab^^- This gives 6 local conditions, so that the two tetrads can now be 
Lorentz-rotated only simultaneously, which leaves only the 6-parameter freedom of local 
rotations in their definition. We then have a simple relation 

1'. = e^ACO^ , (8) 

because 

r.r. = ^>yB^u = e^Ae^^ujB.u;^ = g''" , (9) 
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which reproduces Eq.(jl]). As a result, we can vary the action with respect to and and 
the variation of j^^^ will be obtained by varying Eq.(l8]). In order to take into account the 
conditions ([7]), we add to the Lagrangian the term 

{e'^coB^ - e^a;A^)A^^ (10) 

where A^^ = — A^"^ are the 6 Lagrange multiplies. This allows us to consider and cu^ as 
independent variables. Performing then the variation and expressing A"^^ from the resulting 
equations, we find that A^^ = +A^^. Therefore, A^^ should be at the same time symmetric 
and antisymmetric on-shell, so that it should vanish. 
As a result, we arrive at the following field equations 

rp^-IrSP^= ni^TP + SvrGT^"^^^ , (11) 
K-):n5l = r^m'T^, (12) 



where 



with 



T,^ = r,^-5^A.,, 7r = -^<, (13) 



-/ 



= ill - 3)7^ - 7^7a - I e,,..e-P^'^Ya^K; - j e.^.^e'^^^'Yo^K^^K^K! . (14) 
The Bianchi identities for the left-hand side of Eq. (fTTl) imply the conservation condition 



Vp T{ = , (15) 

(9) 

where Vp is the covariant derivative with respect to gf^u- It is worth noting that the mat- 
ter energy- momentum tensor is conserved independently, in view of the diffeomorphism- 
invariance of the matter action S'(m), 

Vp t'^'^K = . (16) 

{/) 

The Bianchi identities for the left-hand side of Eq.( IT2l) imply that Vp = 0, but these 
conditions in fact follow from Eq. (fT5|) . in view of the diffeomorphism-invariance of the in- 
teraction term ^int = / Cint\/—gd'^x. Indeed, let us consider a diffeomorphism induced by 
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a vector field It induces variations 



(9) (9) 
(/) (/) 

6u^ = -d.CoJ^ - = - V;. Cu^^ - r V. (17) 



while 



= (55int = / <; + 5uj„ \ d X 



(.9) „ , r . (/) 



- / r V. T^V^d'x - U^V. %W-fd'x. (18) 



ig) „ (/) ^„ 

Since can be arbitrary, it follows that the conditions Vp ~ ^ imply that Vp 7^ = 0. 



If — 7- then the source term in f|T2l) vanishes and one obtains 7^^,^ = 0, whose solution 
can be chosen to be fiat metric, which can be written as /^j, = riABdfj,^^du^^ , where 
are sometimes called Stuckelberg fields. All the above considerations then still apply, in 
particular the tetrad formalism, where it is sufficient to choose = (9^$"^. The basic field 
equations are then (fTT|) . (fT5l) . (fT6l) which determine g^^ and 

Let us return to the generic case with 77 7^ 0. In what follows we shall be considering 
solutions of equations (fTT]) . f|T2|) . f|T5|) . f|T6|) with spherical symmetry. We introduce spherical 
coordinates x'^ = (t, r, {}, ip) and choose the tetrads to be 

\ d d did Id Id 

eo=oT^ + ?7r:, ei=p— + -— , 62 = -7^, eg = -— — 77-, (19) 



S dt ^ dr' ' dt N dr' " R d9' ' Rsin^ d^'' 



and 



u° = adt + cdr, u^=ddt + bdr, u'^ = Udi^, = UsmMip, (20) 

where 5", q, N,p, R, a, b, c, d, U are functions of t, r. This implies the spherical symmetry for 
both metrics, while using the residual diffeomorphisms in the t, r subspace one can always 
make the metric g^'^ diagonal. We therefore set 

^°^ = e°eS-eK = |-f = 0, (21) 

so that q = pS/N. Next, we consider the symmetry conditions ([7]), of which the only 
non-trivial one is 

d c 

- e^uif, + e^'wop = e^ul + e^wj = - + qb + ap + — = 0, (22) 
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from where d = —apS — S'^pb/N — Sc/N. We then notice that changing the parameter 
p corresponds to the simultaneous local Lorentz rotations of the two tetrads and does not 
change the metrics. We can therefore impose the gauge condition p = 0, which finally gives 
the following most general expressions for the tetrads: 



eo 



1 A 

Sdi' 



ei 



62 



adt + c dr, 
The corresponding metrics read 



-— dt + bdr, 



]_d_ 
Rde' 

= Ud^. 



63 



1 



d 



Rsin'd dip 
^ - f/ sin -i^dip . 



(23) 



g^^dx^'dx" = S'dt' - N'dr' - R\d^' + sin' M(p 



(24) 



and 



while 



(a2 _ ^) dt" + 2 ^^^^ + *^^^ dtdr - {b' - c') dr^ - U\dd^ + sin' M^^), (25) 



iV2 



N 



7^ 



A« A 



I a/S c/S Q \ 
-cS/N^ b/N 
U/R 
\ U/R ) 

and it is easy to verify that ^^^^'^^ = g^" fau- We also notice that 



-g 



-f 



R^ fab_ 

LP m 



-1 



(26) 



(27) 



We can now compute £int and the tensor rj^ defined by f|T4|) . they are shown in the 
Appendix. Since our fields are only SO (3) -invariant, we have at the time being ^ 0, 
— 7^ 0. Our aim is to find homogeneous and isotropic solutions for (jf^^ , in which case 



and T; - 



one should have T° 
the components r° and — t^- One has 



0. We therefore proceed to eliminate 



R^S 



{-R {3R - 2U) + C3 (3i? -U){R-U)+ a {R - U)^) . 



(28) 



For this to vanish, we can either choose c = 0, or set to zero the expression between the 
parenthesis. 
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III. SOLUTIONS WITH DECOUPLED METRICS 

Let us first consider the case where c 7^ and choose U = CR, where C is a constant. 
Eq.(l28l) then becomes 

rO = I {2C - 3 + C3(C2 - 4C + 3) + c,{C - 1)'}, (29) 

which can be set to zero by adjusting the value of C4, but then one finds 
^ (C - 1)C3 - C + 2 / 2 Ca Cb c2 ab 



which can in turn be set to zero by adjusting C3. It follows that setting 

_C-2 _ C2-3C + 3 

so that €3(03 — 1) + C4 + 1 = 0, one achieves both r ° = and = r^. The r/f components 
shown in the Appendix then reduce to 

= CiC -I) (4, + ^] , (32) 



while 
This gives 
whereas using ([2 



\m NS 



A. = C(C-l)(^ + ^_^). (33) 



T!; = T!:-si:c,,,, = ic-m, (34) 



r/ = -^r,^ = ^5.^ (35) 
The field equations ( |TT|) . ( |T2l) therefore become 

G(; = m2(C-l)5^ + 87rG'T(°')^, (36) 

g^ = r^m'l^Si:^ (37) 

so that the equations for g^,^ decouple from those for /^j,. It is now easy to get cosmological 
solutions. Setting 

S = a(t), Ar = — i2=, i? = ra(t) (38) 

with k = 0, ±1, so that 

g^^dx^dx" = s?{t) {de - - r'(rf^' + sin' Mcp^)^ (39) 
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and choosing SttGT^^^u = diag(p(t), — P(t), —P{t), — P(t)), equations fl36l) reduce to 

3^±^=m\C-l) + p, (40) 

where is defined by the conservation condition 

p + 3-(p + P)=0. (41) 

These equations describes the late time cosmological acceleration. If p = 7P then p ~ 
a-3-3/7 gQ ^]-^g^^ fQj^ large a the second term on the right in f HU]) becomes negligible. The 
dynamic is then driven by the cosmological term ■m'^{C — 1), which we assume to be positive, 
so that C > 1. 

It is worth noting that Eq. fHOj) is exactly the same as Eq.(18) of js] obtained in the RGT 
theory. These solutions therefore do not depend on weather the metric /^i, is dynamical or 
not, which is due to the fact that equations for g^^^, completely decouple from equations 
( l37|l for f^i,. In order to solve equations ( 137|) for f^i, we notice that its components fm = 
and /^^ = f/^sin^i? are already fixed, since U = Cra(t), but /qo, for, frr are still free, 
since they contain three up to now unspecified functions a, b, c. To see that this freedom is 
enough to fulfill the ten equations fl37|) . we notice that one can consider U as the new radial 
coordinate. The time coordinate should also be changed, so that 

t^T{t,r), r^U{t,r), (42) 

and the metric becomes 

f^^dx'^dx'' = fTT dT"^ + 2fTudTdU + fuudU^ - U\d^^ + sin^ Mif'^) , (43) 

where /tt, fru, fuu are functions of T, U. The structure of the source term in ( 1371) does not 
change in new coordinates, so that we should solve the Einstein equations with the nega- 
tive cosmological term rjw?^^^ to find a metric parameterized by the radial Schwarzschild 
coordinate U . The solution is the anti-de Sitter metric 

dip 

f^^dx^dx' = F^dT^-^- {dd^ + sin^ M^^) , (44) 

F 

where F'^iJJ) = 1 + rjm'^^^ U"^. One can now establish the relation to the t,r coordinates, 
since we can read off the tetrad components from (jHj) , but on the other hand they are given 
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by (123|) . so that one can compare to obtain 



This determines 



and also 



w° = FdT = FTdt + FT'dr = adt + cdr 



cu^ = = ^ (adr + rkdt) = — ca/1 — kr^ dt + bdr . (45) 
F F 



6=—, c= , a = FT, (46) 



T 



I FWl - kr^ ' ^^''^ 



Together with U = Cra{t), this estabhshes the correspondence between the t,r and T,U 
coordinates and also specifies all the unknown functions in the solution. 



IV. GENERIC SOLUTIONS 

Let us now return to Eq.(l28ll with arbitrary C3, C4 and set c = 0. This gives r° = 0, while 

r;-r^ = {US-3RS + aR+C3{a-2US + 3RS-2aR)+Ci{-US + Ua + RS-aR)}. 

We now choose both metrics to be homogeneous and isotropic, 

^ = a(t), N= i? = ra(t), a = a{t), b = -M=, U = rf3{t). (48) 

This insures that the energy-momentum tensors depend only on time and have the diagonal 
structure, Tj^ = diag(TO, T;, T;, T;) and %^ = diag{T^, r;,T;,T;) (the explicit form of 
the tensor components can be read off from the formulas given in the Appendix). The 
independent equations are then the two Einstein equations 

G° = mX + P, G'o = , (49) 

as well as the conservation condition for Tj^ , 

T0 + 3-(T0-r;) = 0. (50) 

One can check that the conservation condition for , 

rv 

V + 3-CC-77) = 0, (51) 
a 
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gives exactly the same equation as (l50l) . which shows once again that is identically 
conserved if T^^ is conserved. 

The Gq equation explicitly reads 

k^ + kci^ , ^ , 3/3(3 _ - C4) , 3/32(c4 + 2c3-l) P^Cs + c^) 



- ^4 m ^4c3 + C4-6 + + ]+p 

= m^T^ + p, (52) 

while the conservation condition 

{(3c3 + C4 - 3)a=^ + 2(1 - C4 - 2c3)a/3 + (03 + C4)/3^}(aa - a/3) = 0, (53) 

and the Qq equation 

^/S"^ + ka'^ of 3(c3 + C4)a 3(c4 + 2c3 - l)a2 (3 - 3c3 - C4)a 
^ - rjm [C4 \ — h 



= . (54) 

Let us set the second factor in flS51) to zero, 

a = — , (55) 
a 

thereby solving the conservation condition. Setting /3(t) = cr(t)a(t) Eq. fl52|) reduces to 
a^ -|- /c a^ 

3 = m^{l - or) ((c3 + C4)a^ + (3 - 5c3 - 2c4)cr + 4c3 + C4 - 6)) + p 

Si 



m^T^ + p = p,{a), (56) 



while Eq. flM|) becomes 

Si ~\~ k Si ^ / 9 / \ \ 

3 T = rfm (c4(T — (3c3 + 204)0" + C4 + 3c3 — 3) 

a a 

= r]m^a%° = p,{(j). (57) 
We see that the sources of the two metrics are proportional, 

m^T" + p = r/m Wo^, (58) 
where a = a{p) fulfills the algebraic equation obtained by taking the difference of ( 156|1 and 

(c3 + 04)0-^ + (3 + ?7C4 - 6C3 - 3c4)o-^ + (-9 - 3?7C3 - 37704 + 9c3 + 3c4)o- 

^ r/(3 - 3c3 - C4) ^ _ ^ ^4c3-Q + A,. (59) 
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Since p = p(a) in view of the conservation condition (HTI) . one therefore obtains a = cr(a). 
Injecting this to the right-hand side of Eq.f l57|) (or f l56|) ) gives the source term p*(a), so that 
the solution a{t) can be determined. 

Let us study roots of the quartic equation ( 159|) . first when C3 + C4 7^ 0. For p = 
there are generically two real roots, one of which is cr = 1 with p*(o") = 0, but depending 
on the parameter values there could be altogether four real roots. For example, for 77 = 
1, C3 = —1 and C4 = 4 there are four roots a = —0.93,0.56,1,2.19 with p^{a)/m?' = 
48.78, -0.25, 0, -2.42, respectively. 

For non-zero p there generically remain only two real roots, since the other two, if exist, 
merge to each other and disappear when p increases. When p — )■ 00, one of the two remaining 
roots is defined by 

^KS^Sc^^p^ (60) 

and the second one is 

{c, + c,)a'^A^. (61) 

We shall say that the root (!60l) belongs to the physical branch, since a is small and one can 
see from ( l56l) that Tq = 0(1) and m^lTg | ^ p because m is small, so that p*(p) = p+0{rn?). 
This is physically expected, since the graviton mass contribution to the total energy density 
should normally be small if the matter density is large. At the same time, these natural 
expectations do not apply to the root (|6T]) . since a is then large and 



m^r° = p,-p = -p + 0(p2/3), (62) 

so that the contribution of the graviton mass to the energy is as large as the matter con- 
tribution, and the two actually cancel each other, up to subleading terms. The resulting 
energy density 

pM = + P = C4Vm'a' + 0(a) = , "^"^^/s P"^' + ^ip'^'^ (63) 

|c3 + d^yi-^ 

can even be negative, depending on the sign of C4. We therefore say that the root (ET]) 
belongs to the exotic branch. 

Both the physical and exotic branches p*(p) extend from large to small values of p, so 
that they describe the decrease of p during the universe expansion. When the universe is 
large and p — )■ 0, the total energy p* approaches a constant value that can be positive or 
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negative or zero, depending on tlie parameter values. For tlie pliysical brancli is always 
positive and tends to zero as p —t- if 3 — 3c3 — C4 > 0, while for 3 — Sc^ — C4 < it 
approaches a positive value (for example, p* — t- 15.79 for C3 = —1, C4 = 4, 77 = 1). For the 
exotic branch p* is positive/negative at large p if C4 is positive/negative, respectively, but it 
seems to always approach a non-zero negative value when p — > (if > 0). 









^ — --.^^^^^^^^ ^exotic 





-5 -4 -3 -2 -10 1 

ln(l/p) 




-2 
ln(l/p) 



FIG. 1. The total energy density p* = m^TQ +p versus the matter energy density p for the physical 
and exotic branches for C3 = 0.9, C4 = —1, r\ = \ (left) and C3 = 0.9, C4 = 1, r/ = 1 (right). 



The described above different type behaviour of p*(p) can be seen by solving the algebraic 
equation (!59l) numerically for different parameter values, as shown in Fig.l. One more type 
of solutions, shown in Fig. 2 (left panel), is obtained by changing the sign of r], in which 
case the sign of p* along the exotic branch changes from negative to positive values as p 
decreases. 

If C3 + C4 = then the coefficient in front of the highest power in Eq. fl59|) vanishes, so 
that there remain three roots. One finds in this case three different branches p*(p), these 
are the physical branch ( l60l) and two exotic branches that start at large p when one replaces 

m by 

(3 + r/C4 + 3c4)a2 = , (64) 

since there are two possibilities to choose the sign of a when one takes the square root (these 
branches are called in Fig. 2 exotic+ and exotic—). If p is large, then the energy for both 
exotic branches is the same up to subleading terms. 



p*(p) = m^Tg° + p = Cij]m^a^ + 0(cr) 



7704 



.p + 0(//^). 



3 + 7704 + 3c4 

The behaviour of p*(p) in the whole range of p is shown in Fig. 2 (right panel). 



(65) 
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^exotic— 



-2 

ln{l/p) 



FIG. 2. The total energy density versus p for the physical and exotic branches for C3 = C4 = 1, 
r/ = — 1 (left) and C4 = — C3 = 1, r/ = 1 (right). 

Now that we have determined p*(p), we can proceed to find solutions for the scale factor 
a(t). Introducing the physical time dr = SL{t)dt, the Einstein equation fl57p becomes 

which describes a 'particle' with the total energy —k = 0, ±1 moving in the potential 



U(a) = -y p*(p). 



(67) 



Assuming a specific equation of state P = P{p) for the matter, the conservation condition 
(HTj) gives p(a), which allows us to compute U(a). In Fig. 3 we show U(a) computed with 
the ultra-relativistic equation of state, 

p = 3P ^ P(a) = ^. (68) 

Solutions of Eq.f l66p for a given k correspond to the regions of a where U(a) < —k. Inspecting 
the U(a) curves in Fig. 3 reveals then five different cases, of which the first two (in the left 
panel) correspond to the physical branches. 

Type I solution correspond to the physical branches with p*(p) — ?■ for p — )■ 0, in which 
case U(a) is negative and tends to zero as a — )■ 00. The scale factor a(r) behaves qualitatively 
in the same way as in the matter dominated universe: it ranges in the finite limits in the 
spatially closed case = 1, it linearly grows with r for k = —1, and it increases as y/r for 
A; = 0. 

Type II solutions correspond to the physical branches with p*(p) — )■ p*(0) > (as for 
example in the right part of Fig.l). At early times they coincide with the ordinary matter- 
dominated cosmologies, since for the physical branches one has p^, ~ p if p is large. However, 
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for large a one has U(a) = — p*(0) and da/dr ~ a, so that at late times solutions with 
k = 0,-1 enter the phase of accelerated expansion. For k = 1 the things are slightly more 
subtle. The three curves Ila, lib and lie in Fig. 3 correspond to different choices of the 
integration constant po in (EH])- If Po is large (curve lie) then U(a) < —1 and the solutions 
are similar to those with k = 0, —1. If po is small then the potential can exceed the value — 1 
(curve Ila). Then there is a solution for which the 'particle' rests on the left of the reflection 
point A (see Fig. 3) close to the cosmological singularity a = 0, but there is also the solution 
that stays on the right of the reflection point B (see Fig.3), it never approaches singularity 
and shows acceleration at large a. 



1 

0.5 
g 
Is" -0.5 
^ -1 
-1.5 
-2 
-2.5 
-3 





k=-1 








k=0 










\IIa 






k=1 / 






/ ''^ 














A lie 


0.2 


0.4 


0.6 


0.8 



2 
1.5 
1 

0.5 


-0.5 
-1 
-1.5 



1 iv\ 


/ill ' / 


/ill 


^^11- 
















k=-1 








k=0 








k=1 





1 

a 



FIG. 3. The effective potential U(a) ()67p for the physical (left) and exotic (right) branches for 
solutions shown in Fig.l and Fig. 2. 



The exotic solutions are shown in the right part of Fig.3. Type III corresponds to exotic 
branches for which p*(p) is positive at large p but becomes negative when p is small. The 
potential U(a) then grows from minus to plus infinity and the 'particle' is always confined 
to the region close to singularity. Type IV corresponds to the exotic branches for which 
p^, is always negative, so that U(a) is unbounded from above but has a positive minimum. 
Solutions can exist only for k = —1 and describe oscillations in the potential well between 
the two reflection points {C,D in Fig.3). Finally, type V corresponds to the exotic branch 
for ?7 < for which p^, changes from negative to positive values as p decreases. The potential 
U(a) is then monotonically decreasing (see Fig.3) and the solutions always stay away from 
singularity and show the self-acceleration at large a. 

Summarizing, only types II and V show self-accelerating solutions. Type II solutions 
arise in theories with C4 > and C3 + C4 7^ 0, they evolve as the matter-dominated universe 
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at early times, but enter the accelerated phase at late times. Type V solutions show late 
time acceleration, while at early times they are regular, being repelled from the singularity 
by the negative total energy p*. We notice, however, that such solutions require the second 
gravitational coupling constant to be negative. 



V. NON-ACCELERATING SOLUTIONS 



Let us now return to the conservation equation ( l53l) and try to fulfill it by setting to zero 
its first factor and not the second one. We therefore abandon the condition ( |55l) . but require 
instead that /3(t) = aa{t) where cr is a constant. Then Eq. (l53|) will be satisfied if 



(c3 + C4)a2 + 2(1 - C4 - 2cs)a + Scs + C4 - 3 = 0, (69) 



Eq. (!52!) then reduces to 



a4 



(4c3 + C4 - 6 + 3a(3 - 3c3 - C4) + 3a'^{c4 + 2c^ - 1) - a^{c^ + C4)) + p 
A(a)+p, (70) 



while Eq. (l5^ becomes 

f _ 3(03 + 04) ^ 3(C4 + 2C3 - 1) ^ (3 - 3C3 -c,) . 



(71) 



Combining (ITOjl and (ITTj) one obtains 



2 2 2(A(^) + p)a2-3A; 



rjm?a'^a?TQ —3k ^'^ ^ 

The quadratic equation fl69|) has two roots, 

_ 203 + 04 - 1 lb a/o3(o3 - 1) + 04 + 1" , , 

C3 + C4 

and the value of the cosmological constant A(o") in Eq. (!70|) is positive for one of them and 
negative for the other one. There is, however, an additional condition, since must be 
positive, and as the numerator in (172]) is positive because (A((t) + p)a^ — 3k = 3(a/a)^ > 0, 
it follows that Tq should be positive, since otherwise the denominator will become negative 
for large a. Now, it turns out that if 7^° > then A (a) < 0, while if A(cr) > then 7^° < 0. 
This eliminates solutions with A(cr) > 0, but there remain solutions with A(cr) < and 
%° > 0, which exist if A; = 0, -1. 
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VI. LIMIT r] 







Let us first consider the solutions with decoupled metrics of Sec JIIII Taking the limit 
?7 — )■ does not affect the physical metric g^i, determined by Eqs.(l39l). ( I40l) . ( 14T1) . On the 
other hand, the metric ff^u in (jHj) becomes flat, since F — t- 1 when r/ — )■ 0. One can write 
f^,u = VABdf.^'^d^^^ with <l>° = T{t,r) and = U{t,r)n'' where f/(t,r) = Cra(t) and 
= (sini^cos*/?, sint^sinyj, COS"*?), while T{t,r) is obtained by settin F = 1 in (j4j 
Cr2 



T(t,r) 



a if k = 0; T{t,r) = kCaVl - kr^ if k = ±1. 



(74) 



Equation ([39]) , (gO]) , gl]) , (EH) exactly agree with Eqs.(16)-(18) obtained in [sl in the RGT 
limit. We therefore conclude that the solutions with decoupled metrics have the counterparts 
in the RGT theory, to which they approach when — )■ 0. 

Let us now consider the generic solutions of SeclIVI A direct inspection shows that /^i, 
does not necessarily become fiat when 77 — 0, because the source term rjm'^TQ in Eq. (|5^ 
does not then vanish, neither does the source for the physical metric = r^m^cr^T^^. To 
understand how this is possible, we notice that for the physical branches a is small when p is 
large, because rj/cr p/rn? (see Eq.( l60|) ). On the other hand, Eq.( l57|) shows that when cr is 
small then p^, ~ rjim? /a = p + . . . for any rj. As a result, the effective potential U = — a^p*/3 
does not vanish in the region where a is small but approaches a non-trivial limit as 77 — )■ 
(see Fig. 4). 



OS 





FIG. 4. The effective potential U(a) (j67p for the physical solutions with C3 = 0.9, C4 = 1 (left) and 
exotic solutions with C3 = 0.9, C4 = — 1 (right) for several values of rj. 



For some exotic branches a becomes very small for small t] when p — i- 0, in which case 
one finds from (|59|1 a ~ ri{3 — 803 — C4)/(c4 + 4c3 — 6) so that p^ ^ rj/a is independent of 
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7] and the effective potential U{a) does not vanish at large a. For other exotic branches a 
never approaches zero, in which case — > as r/ — > 0, therefore both g^^i, and f^^, become 
flat. The conclusion is that the generic solutions of Sec JIVI do not have non-trivial analogs 
in the RGT limit. 

Finally, for the solutions of Sec JV] the metric g^i, is determined by ( 169|) . ( 1701) and does not 
depend on rj. The source term in (|7T1) vanishes for — t- and f^,, becomes flat. The limit is 
possible only for k = —1, since a in ( 1721) becomes ill-defined for A; = if — i- 0. As a result, 
the solutions of SecJV]do have, for k = —1, analogs in the RTG limit. Moreover, for r] = 
one can choose A(cr) > in (!70|) . as this no longer contradicts the positivity of in (172|) . 
Such solutions were found in 9|, but only for A(cr) < they can be extended to r] ^ 0. 

We see that the bimetric theory admits solutions which do not approach for r/ — )■ those 
of the 1] = theory, and vice versa, the 1] = theory has solutions which do not generalize 
for 77 7^ 0. One can construct more solutions for 77 = if we go directly to Eqs.(l52l).( l53i) and 
require that the metric parameterized by the functions a, /3 (with U = rP{t)), 

1 — kv^ 

be fiat. It will be fiat if one finds T{t,r) such that 

dT^ - dU^ = a{tfde - dr^ , (76) 

1 — kr^ 

which is equivalent to three conditions 

^2_y,2^^(^^ ^2 _ ^2^2 ^^2^ fr = r$/3. (77) 
1 — kr'^ 

One possibility to fulfill these conditions is to set a = 0, /3 = C, and 

C 



0: T = 0, k = ±l: T = ^= Vl - kr^ . (78) 



The conservation condition (135]) is then fulfilled and one is left with the Einstein equation 
(!52|) where /3 is constant. This reproduces the solutions given by Eqs.(19),(20) in |8| (the 



opposite sign convention for C3 is used in jsj), they exist only in the RGT limit and do not 
generalize for t] ^ 0. 

Another possibility to fulfill ( I77|) is to choose k = —1 and set 

T = VTTr^/3{t), a = p. (79) 
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With f3{t) = aa(t) the conservation condition (153|) is fulfilled if a is given by (173|) . while a 
is then determined by f lTOj) . The solutions were obtained in Ref. jol, they generalize to rj ^ 
if only one chooses the root of (1731) for which A(cr) in ( ITOl) is negative. 

Summarizing, among the accelerating solutions of the RGT theory only the special so- 
lutions (|74l) generalize for r] ^ 0, while among accelerating solutions of the bimetric theory 
only solutions with the decoupled metrics of Sec JIIII have the RGT limit. 

Recently it was claimed in the literature that the RTG theory does not actually admit 
homogeneous and isotropic cosmological solutions 12| (apart from those obtained in the 
decoupling limit 13j). At the same time, the presented above analysis shows very explicitly 
that such solutions exist, thus confirming the results of 8|, {q]. The negative argument of 12 | 
assumes that in the unitary gauge, where = x'^ and f^^, = t]^^, the physical metric g^i, is 
diagonal (see Eq.(13) in [12]). However, the two metrics cannot in general be diagonal at the 
same time. For example, g^^, is diagonal in coordinates t, r, {}, ip, but f^^, = rjAsd^^^dy^^ 
with $^ defined by formulas around Eq.f l74p is not diagonal. For the solutions (1781) both 
metrics are diagonal at the same time, but f^u is degenerate, so that the argument of |l2 | 
again does not apply. 



VII. SUMMARY 

We have presented the homogeneous and isotropic cosmological solutions within the bi- 
metric generalization of the RGT massive gravity theory. These solutions can be spatially 
open, closed, or fiat, and at early times they are sourced by the perfect fiuid, while the 
graviton mass typically manifests itself at late times by giving rise to a cosmological term 
whose value is determined by the theory parameters 03,04,77. In addition, there are also 
exotic solutions for which already at early times, when the matter density p is high, the con- 
tribution of the graviton mass to the energy density is large and screens that of the matter 
contribution. The total energy m^Tg + p can be negative, which can lead to non-singular 
solutions, as in the case of type II solutions with k = 1 oi Sec JIVI For type V solutions 
of Sec JIVI the cosmological singularity is removed altogether, but this requires the second 
gravitational coupling to be negative. 

In the limit where the second gravitational coupling tends to zero the generic solutions of 
Sec JIVI do not reduce to solutions of the RGT theory, since both metrics remain then curved. 
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However, the special solutions with decoupled metrics do have the non-trivial RGT limit. 

The analysis of stability of our solutions remains an open issue to study. Since the graviton 
contribution to the total energy can be negative and very large for the exotic solutions, it 
is not impossible that the ghost could be still present in the theory, which may affect the 
stability. 



APPENDIX. 



Here we list the energy momentum tensor components in the spherically symmetric case. 
Using the expression ( 126|) for 7^^^ and computing Kf} = 5^ — 7^^^ gives the following value of 
the interaction Lagrangian ([2]): 

ah 3a 36 2aU 2bU QU 

R-U f2ab 3b 3a Ua 2U bU 



R \NS N S RS R RN 

iR-Uf f a b ab ^\ , , , 

while the non-zero components of defined by Eq. f|T4|) read 

ab 3a 2aU R-U (3a 2ab 2c^ aU 



SN S SR R \S SN SR 

-Uf (a ab 
~R2 \'S ^ 



(R-Uf fa ab , 



ab 3b 2bU R-U (3b 2ab 2c^ bU \ 

^ R^ \N SN mj' 



Ufa b _ U 

r[s^n'^^r 



+ C4 



U{R - U) 
^2 



U ( 2b 2U bU 



aU ab 



2a 

~s^Jr^ 



a b 



ab 



(A.4) 



{-R {3R - 2U) + C3 {3R -U){R-U)+ a {R - Uf) . 



(A.5) 
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The components of the two energy- momentum tensor are then simply obtained from Eq.( !T3|) . 
where ^J^^ / y/—f is given by fl27|) . 
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